Gravitational field measurement with an equilibrium ensemble of cold àtoms 
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A new approach to the measurement of gravitational fields with an equilibrium ensemble of 
ultra-cold alkali àtoms confined in a celi of volume V is investigated. The proposed model of the 
gravitational sensor is based on a variation of the density profile of the ensemble due to changing 
of the gravitational field. For measurement the atòmic density variations of the ensemble the 
electromagnetically induced transparency method is used. 
PACS number(s): 39.20.+q, 03.75.Dg, 04.80.-y, 32.80.Pj, 42.50.Gy 



Gravitational sensors operating on tcchnologies of 
atom òptics usually use a variety of atom interferometry 
techniques for the detection of gradient of the gravita- 
tional potential Atom interferometry requires ele- 
ments equivalent to beam splitters and mirrors of light in- 
terferometry • These elements complicate the grav- 
itational sensor's design and impose limitations on its 
sensitivity. The objective of this work is to provide a 
scheme of gravitational sensor excluding atom interfer- 
ometry technique. We will investigate the posibility of 
gravitational field measurement with an equilibrium en- 
semble of ultra-cold alkali àtoms (medium) confined in a 
celi of a volume V. The density profile of this medium 
is changed due to the interaction of its àtoms with the 
gravitational field. The gravitational field can be mea- 
sured by measuring the phase shift of an optical beam 
passing through the medium. For measuring the density 
variation of the medium, the electromagnetically induced 
transparency technique is used p 7 ]- 9J . The sensitivity of 
the idialized model of the proposed gravitational sensor 
is comparable with the sensitivity of modern atom inter- 
ferometric gravimeter. 

Let us consider an atòmic cloud confined and cooled 
inside a thin-walled celi of a small volume V. The tech- 
nique of alkali atom cooling is well developed and pre- 
sented in a number of publications 0|"E3· After the 
trapping potential is turned off the cloud will expand 
filling all volume V. We will consider that on this stage 
the atòmic ensemble is only under the influence of a grav- 
itational field and will neglect the heat exchange between 
the ensemble of alkali àtoms and the celi during the mea- 
surement time. The field theoretical Hamiltonian of the 
ensemble in the presence of the external gravitational 
field is written as 



-— V 2 + V'(f) 
Im 



ip(x, t) 



where tp(x, t) and tp + (x, t) are field operators, V(x) is the 
atom-gravitational field interaction potential, m is the 



atòmic mass, and U (x, x') is the two body interatomic 
potential [12|. Operator ift(x,t) creates an atom at the 
point x and time t while tp + (x, í) annihilates an atom. 
The field operators obey the following commutation re- 
lations: [ip(x, t), ip(x', i)] = and [4>{x, t), ip + {x\ t)] = 
5 3 (x — x 1 ), which define the statistics of the ensemble. 
We can also introduce a number operator 



N = d 3 xip + (x,t)ip(x,t) 



(2) 



with the eigenvalues interpreted as the number of the 
àtoms. For further analysis we will use the grand parti- 
tion function of the atòmic medium 



Z = Sp exp (0(pN - H)) 



(3) 



where fj, is the chemical potential, and f3 = 1/ksT is 
the inverse temperature {ks is the Boltzmann constant). 
If àtoms of the medium are under the influence of only 
a gravitational field, then the potential of interaction is 
given by the classical equation V(x) — m<&{x), where 
<£>(;?) is the Newtonian gravitational potential. Inserting 
Q and © into © and neglecting interatomic interac- 
tions |l3| we obtain the grand partition function in the 
form 



Z = [][!- exp (-/3(e,- /i))]" 



(4) 



where q is the energy of a boson in z-th state. Statisti- 
cal and thermodynamical properties of the medium can 
be investigated with the help of a grand thermodynamic 
potential given by O = — (l//3)lnZ. Using the equation 
j|4*|l we obtain the expression for the thermodynamic po- 
tential in the form 



(1//3) ]T ln [1- exp (-/3( ei -/j))] 



(5) 



In quasi-classical regime we can replace in the equation 
|J5J| the summation over single atòmic states by the inte- 
gration J2i ~> 3 J (2nh)~ 3 dT , where dT = d 3 pd 3 x, and s 
is the number of the possible spin states of the atom. The 
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atòmic energy in this case can be represented by a sum 
of kinetic and potential energies = p 2 /2m + m$(x) 
and the expression for the grand thermodynamic poten- 
tial will be 



fi = (V/3) 



sdT 



ln [l - z($) exp(~(3p 2 /2r 



(6) 



where the gravity dependent fugacity z is given by the 
equation 



z($) = exp [/3(p — m$)] 



(7) 



The grand thermodynamic potential (jüj) can be written 
in the form íï = J d 3 xuj(x) where the density of grand 
thermodynamic potential is given by the equation 



w = (1//3) 



sd 3 p 



ln [l - z($) exp (-/3p 2 / 2 ™)] ■ (8) 



The average density of àtoms in the medium is found as 
N = —(dui/d/j.) and has the following form 



TV = 



(9) 



(2tt^) 3 exp(/3p 2 /2m)) - l' 

The integration of JÜJ yields 

^ B = 53/2 (z), (10) 



where XdB — y2Trti 2 (3/m is the thermal de Broglie wave- 
length and 53/2(2) is a function defined by the equation 



9k{z) = 



n=l 



The functions gk{z) obey the recursion relations 



z -r9k{z) = gk-i(z) 
dz 

and have the following integral representations 



(11) 



(12) 



9k{z) 



r(fe-i/2) 



V^F r(2fc- 1) y 



„2fc-l 



exp(x 2 ) — 1 



(13) 

where T(fc) is the gamma function. The functions gk{z) 
have singularity at z = 1 , which leads to BEC, moreover, 
the series llllfl converges only for k > 1. The equation 
(|10fl describes the relation between the density of parti- 
des N, temperature T and the gravity dependent fugac- 
ity z. The variation of the atòmic density of the medium 
is found from Í|1C)|> in the form 



SN = -f3mN F(z )S<S>, 



(14) 



where F{z) is given by the equation F(z) = 
9i/2{z) I '33/2 (-z)- I n the equation (|14|l the atòmic density 
variation SN = N—No, an d the variation of gravitational 




1 driving beam 

FIG. 1: Mach-Zehnder interferometer deteets the phase shift 
of the probe beams after they pass through the cooled atòmic 
medium placed in the celi of the volume V — 1 cm 3 . The 
gradient of gravitational potential causing the phase shift is 
given by an arrow. The celi has the length L, and the upper 
and the lower beams of the interferometer are separated by 
the distance Y. Dashed arrows correspond to driving fields. 



potential <5$ = $(a?) — $(5?o) are defined between points 
x and Xq. Fugacity zq at point xq is related to iVo and 
temperature T of the medium by the equation (|10fl . Lct 
us select the coordinate system {o; 1 , x 2 , x 3 }, where the 
vector of gravitational acceleration is g = {0,0, —g}. In 
this coordinate system the variation of gravitational po- 
tential can be written as SQ = gY, where Y = (x 3 — Xq) 
is the displacement between two points along x 3 -axis. 

The equation (|14f) can be used for the analysis of the 
influcnce of gravitational field on the medium consisting 
of non-interacting alkali àtoms in the celi of volume V. 
Lct us rewrite it in the form 



SN = -N f3mYF(z )g- 



(15) 



As follows from this equation, the gradient of gravita- 
tional potential causes the variation of the atòmic density 
of the medium in the celi that can be detected with elec- 
tromagnetically induced transpareney (EIT) technique in 
the combination with the method of optical interferome- 
try (Fig. 1). 

Let us consider EIT effect in application to the medium 
in gravitational field by taking into account the infernal 
structure of bosons. We will assume that each boson of 
the medium is a three-level system with one upper (a) 
and two lower levels (b) and (c). To obtain EIT con- 
ditions, two optical (driving and probe) fields are used. 
Let the driving field of Rabi frequeney íl couple the levels 
(a-c) whereas the probe field of frequeney u couple the 
levels (a-b). The driving field provides an interference of 
different possible absorption pathways producing no ab- 
sorption effect 14]. As the result, the medium becomes 
transparent to a probe field and the real part of suscep- 
tibility x' of the probe frequeney near the resonant point 
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will be 



X 



he ft 2 ' 



(16) 



where p is the o — > b transition matrix element, N is the 
atòmic density, and A a {> = uj^ — lo is the atom-probe field 
detuning [15J . The dispersive properties of the medium 
are dehned by refractive index n = 1 + x'/% which is 
writtcn in the form 



n=l- 



A 



2he n 2 



N. 



(17) 



As follows from Fig.l, the interferometric phase shift <fi 
of probe beams is a function of the variation of refractive 
index Sn of the medium 



2nL 



Sn, 



(18) 



where A is the wavelength of the probe field, and L is the 
length of the celi. Interaction of àtoms of the medium 
with gravitational field changes the density profile of the 
medium. The density variation in the medium leads to 
variation of the refractive index along the gradient of 
gravitational field. When two separated probe beams 
pass the medium at points of different gravitational po- 
tential, they gain the phase shift depending on the gradi- 
ent of the gravitational potential. Therefore the gradient 
of the gravitational potential can be measured by measur- 
ing the phase shift of probe beams in the interferometer. 

From the equations i|15|) . (|17fl and l|18l) we obtain the 
relation between the phase shift and the gravitational ac- 
celeration. The variation of the refractive index is found 
from the equation (|17|) as 



Sn 



3 A 3 7 A 
lir 2 n 2 



SN, 



(19) 



where 7 = p 2 uj 3 / (6irheoc 3 ) is the radiative decay from 
the upper lcvcl. Thcn the phase shift of the two beams 
in the interferometer will be 



3 L\ 2 -fA 

4^ n 2 



SN. 



(20) 



Inserting the equation l|15|l into the equation l|20|l we get 

3 A\ 2 j 2 xN Q mF(z )~ 



4tt 



k R TVL 2 



9, 



(21) 



where A = LY is the area of the side surface of the 
medium, and x — A/7 is the dimensionless detuning 
seale. The limiting sensitivity of the device is dehned 
by the equation 



3 AX 2 j 2 xN mF(z ) 



4tt 



k B TÜ 2 



(22) 



For numerical demonstration of the limiting sensitivity 
of this gravitational sensor, we will assume that the driv- 
ing held is strong and Í2 « 7 |ïij. Taking into account 
that F(z) > 1, and putting for simplicity xF(z) w 1 we 
can rewrite the equation (|22|l in a simple form [17| 



Sg = 



4tt 



3 AX 2 N m 



TS<t>. 



(23) 



For a wavelength A = 500 nm, atòmic density Nq — 
10 12 cm~ 3 , atòmic mass m — 10 -26 kg, and the side 
area of the celi A = 1 cm 2 , we fmd that the coefhcient 
of the equation (|23|l in square bracket is of the order of 
ten and, therefore, one can simplify the equation (23) 
as Sg w lQTScf). Assuming the temperature of the gas 
is about T = 1 mK and computing the sensitivity of 
the interferometer with the equation S(f) — yjhui/ P p t for 
P p = 1 mW probe beam power and r = 1 second of 
measurement time, we fmd that the minimal detectable 
gravitational held changing is Sg sa 10~ n g. As follows 
from the equation i|23|) , the important factor contributing 
to the sensitivity of the sensor is the side area of the celi 
A. The volume V of the celi contributes to the sensitiv- 
ity due to changing the density of àtoms No- Therefore 
probe beam separation Y (height of the celi), length of 
the celi L and its thickness can be adjusted based on 
optical interferometer design. 

In conclusion, the numerical analysis has demonstrated 
that cold-atom EIT gravitational sensor operating with 
an equilibrium ensemble of cold alkali àtoms confined in 
one cúbic centimeter volume allows us to measure gravi- 
tational acceleration with high sensitivity. The sensitiv- 
ity of the proposed gravitational sensor operating with 
mK temperature alkali medium is about 10 _11 g which is 
higher for the same operational time then that of modern 
cold-atom gravimeters operating on atom interferometry 
effect. Lowering the temperature of alkali medium to //K 
regime, which is the operational temperature of the atom- 
cloud gravimeter, one can reduce the number of àtoms 
from No — 10 12 to Nq — 10 9 cm~ 3 keeping the sensor's 
sensitivity 10~ n g. As follows from the above analysis, 
obtaining high sensitivity requires high atòmic density 
and precise control of cloud temperature. The better 
we control temperature the longer is the time of phase 
measurement. The idialized model which has been inves- 
tigated here does not include thermal exchange between 
alkali cloud and the celi. We assumed the temperature is 
constant during the measurement time r. Thermal inter- 
action shortens the time of measurement and lowers the 
sensitivity of the sensor. A possible way to reduce the 
effect of thermal exchange on the sensor's sensitivity is 
to compensate for the raise of the temperature by experi- 
mentally dehned (or modelled) function T(t), t G r in the 
equation Ij23(l . Due to time dependence of the tempera- 
ture of the ultra-cold cloud the proposed model operates 
in discrete regime. For each conseqüent measurement the 
celi must be rehlled with a new collection of àtoms. 
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